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We show that Parity- Violating Deep Inelastic Scattering (PVDIS) of longitudi- 
nally polarized electrons from deuterium can in principle be a relatively clean probe 
of higher twist quark-quark correlations beyond the parton model. As first observed 
by Bjorken and Wolfenstein, the dominant contribution to the electron polarization 
asymmetry, proportional to the axial vector electron coupling, receives corrections 
at twist-four from the matrix element of a single four-quark operator. We reformu- 
late the Bjorken/ Wolfenstein argument in a matter suitable for the interpretation of 
experiments planned at the Thomas Jefferson National Accelerator Facility (JLab). 
In particular, we observe that because the contribution of the relevant twist four 
operator satisfies the Callan-Gross relation, the ratio of parity- violating longitudinal 
and transverse cross sections, TV Z ', is identical to that for purely electromagnetic 
scattering, TV ', up to perturbative and power suppressed contributions. This re- 
sult simplifies the interpretation of the asymmetry in terms of other possible novel 
hadronic and electroweak contributions. We use the results of MIT Bag Model cal- 
culations to estimate contributions of the relevant twist four operator to the leading 
term in the asymmetry as a function of Bjorken x and Q 2 . We compare these esti- 
mates with possible leading twist corrections due to violation of charge symmetry in 
the parton distribution functions. 
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I. INTRODUCTION 

In the 1970s parity-violating deep inelastic scattering (PVDIS) of longitudinally polarized 
electrons from deuterium played an important role in confirming the Standard Model (SM) 
of particle physics [IH3] . The asymmetry 

A RL = 2LZ£L, (1) 

or + °l 

with a^L corresponding to the scattering cross-section with positive and negative helicity 
electrons respectively, is an excellent probe of the parity-violating electroweak interactions 
in the SM. The results of measuring this asymmetry in the early experiments at SLAC 
led to the correct description of neutral weak interactions well before the discovery of the 
Z boson at CERN and provided a measurement of the weak mixing angle sin 2 9w Since 
then parity- violating electron scattering from various targets has been studied at JLab [IH7] , 
MIT/Bates [SJ [H], Mainz [TTJJ E], and SLAC [12] as a tool for probing physics beyond the 
SM and hadronic structure. 

Currently, an active program is underway at JLab to continue these studies with a new 
level of precision. The Q-Weak experiment [ 1 3 J , which will measure the weak charge of 
the proton at low electron momentum transfer (Q 2 ), is expected to determine sin 2 9w to 
0.3% precision, making it the most precise test of the running of the weak mixing angle to 
date. Furthermore, the recently approved 12 GeV upgrade of CEBAF at JLab, expected 
to be completed by 2014, aims to begin the next generation Moller and electron-deuteron 
scattering experiments. The SOLID proposal [Hj for precision parity-violating electron- 
deuteron scattering, approved as part of the 12 GeV upgrade, will measure A RL over a wide 
kinematic range in Q 2 and Bjorken-x to within 1% at each kinematic point. In addition, 
one high-precision PVDIS experiment with deuterium has completed data taking at selected 
kinematic points with the 6 GeV [T3] beam and another is approved to run at the 12 GeV PS] 
beam. These present and prospective high-precision experimental measurements present new 
challenges for their theoretical interpretation. In particular, substantial uncertainties in the 
theoretical interpretation of the deep inelastic asymmetries will remain unless various effects 
contributing to the asymmetry such as new physics beyond the SM, sea quark distributions, 
Charge Symmetry Violation (CSV), and higher twist contributions are well understood and 
disentangled from each other. Addressing one aspect of these issues is the subject of this 
paper. 

The theoretical interpretation of the deuterium asymmetry can be facilitated by express- 
ing it in the following form 



A 



RL 



( G F Q 2 \ 
\ Ay/2na j 



f\ , 9Vv *i 



F? ' 2 ~ 3 F[ 
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(2) 



Here, g v (g^) are the vector (axial vector) couplings of the Z-boson to the electron; , F? z , 
and F^ z are the structure functions arising, respectively, from hadronic matrix elements of 
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the vector electromagnetic (EM) current, interference of the vector EM and vector weak 
neutral current (WNC), and interference of the vector EM current and axial vector WNC; 
and Y 13 are functions of kinematic variables and the ratios R 1 and W z of longitudinal 
and transverse cross sections for purely EM and WNC-EM vector current interference cross 
sections. In the SM, at leading twist and in the absence of CSV effects, the Y 1 term in Eq.Q 
is independent of y and depends only on g\ and the vector current coupling of the Z-boson 
to quarks [3]. Since g v = — 1 + 4sin 2 6V ~ —0.1, the Yi-term dominates the asymmetry, 
making its scrutiny particularly important for the interpretation of the Jefferson Lab PVDIS 
program. 

Considerable theoretical effort has been devoted to disentangling the various contribu- 
tions to the asymmetry. The effect of twist-four contributions to the asymmetry was first 
considered in papers by Bjorken and Wolfenstein [XT] ITS] more than thirty years ago, where 
it was shown to arise in the dominant axial electron coupling term from a single, non-local 
four-quark operator in the limit of good isospin, negligible sea-quark and CSV effects, and 
up to corrections in a s (Q 2 ). Quantitative estimates of twist-four effects were first obtained 
in [19] where the contribution of the spin-two operators was estimated using the MIT Bag 
Model. This analysis was extended in [20J to include corrections to the F^ z structure func- 
tion. More recently, higher twist effects to the asymmetry were estimated by the authors 
of Ref. |21] , who considered the possibility that W and RJ Z receive substantially differ- 
ent contributions from finite-Q 2 effects. These authors argued that such a difference could 
introduce hadronic uncertainties that might impede the extraction of CSV effects from Am,. 

In this paper, we draw on the observations of [17, 18J that the twist-four contribution to 
the Y\ term in Arl for deuterium, given in Eq. (pj), arises from a single four-quark operator 
involving up- and down-quark fields 

GEO*) = l[u(x)Yu(x)d(0)rd(0) + (it d)} (3) 

to revisit the analysis of Ref. (21]. Noting that the contribution of 0^ d (x) to the electroweak 
structure functions satisfies the Callan-Gross relation at leading order in the strong coupling, 
we find that 

W z = W and Y 1 = l, (4) 

at twist-four up to perturbative corrections. Consequently, all twist-four effects entering the 
dominant term in the asymmetry reside in the ratio F^ Z / 'F± . 

Using the power law dependence in Q 2 of the twist-four effects to the l^-term it may be 
possible, with the precision and the wide kinematic range of the PVDIS program at JLab 
and its possible extension at an electron-ion collider [22J , to disentangle twist-four effects 
from CSV effects depending on their relative overall sizes. To provide theoretical guidance 
for such a program, we utilize the MIT Bag Model [23] to estimate the size and variation of 
the twist-four contribution with Bjorken-x and Q 2 . These estimates extend the earlier work 
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of Ref. [20] by allowing for the ^-dependences of the twist-two and twist-four contributions 
to F^ 1 to differ. We find that if the MIT Bag Model reasonably estimates the magnitude 
of the twist-four contribution from 0^(x), the impact on the asymmetry would likely be 
too small to be extracted without further improvements in experimental precision. In this 
case, however, the planned PVDIS experiments could in principle provide a theoretically 
clean probe of possible contributions from CSV and/or physics beyond the SM. Conversely, 
the observation of significant power corrections to the Y\ term would signal the presence 
of relatively large and theoretically interesting quark-quark correlation contributions to the 
electroweak structure functions. 

Our analysis leading to these conclusions is organized in the remainder of the paper as 
follows. In Section [TTJ we provide an overview of the structure of the deuterium asymmetry, 



setting the context for our analysis of the twist-four contributions in Section III In Section 



IV] we give our MIT Bag Model estimates and compare these with recent parameterizations of 
CSV contributions as well as possible effects from "new physics" in Section [Vj We summarize 
our conclusions in Section VI In appendix |Aj we also recast the argument of [24TT26] in the 



language of the Soft-Collinear Effective Theory(SCET) [2TH29] that shows manifestly that 
the twist-four matrix element contributing to the l^-term satisfies the Callan-Gross relation 
at tree level in the matching. 

II. OVERVIEW 

Before presenting the formalism and derivation of our results, we provide an overview of 
the structure of the deuterium asymmetry and the context for the higher twist contribu- 
tions. The SM parity violating interactions of the electron with the quarks, obtained after 
integrating out the Z-boson, are parameterized as 
G F 



C 

y/2 



e^^ 5 e(C lu u^^u + C ld d^^d) + eYe(C 2u u^^ 5 u + C 2d d^ ll 'y 5 d) , (5) 

where the coefficients C\ q and C\ q are given by 

C lq = 2p NC It (/f - 2Q q k sin 2 6w) (6) 

C 2q = 2p NC II (l 3 e - 2Q e £sin 2 Ow) - l -\ q 2 . (7) 

Here is the third component of weak isospin for fermion /, Qj is the electromagnetic 
charge, and 9w is the weak mixing in the MS scheme. The quantities Pnc, and encode 
the effects of electroweak radiative corrections and take on the values one, one, and zero, 
respectively, at tree-level, leading to 

CT = - \ + \ sin 2 6 W , CtT = \ - \ sin 2 9 W , 

Cf: e = - \ + 2 sin 2 9w, ClT = \-2 sin 2 9 W . (8) 
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The reason for the high sensitivity of A^l to these interactions is that in the limit of 
good isospin and negligible sea quark effects, all hadronic effects are known to cancel in the 
asymmetry at leading order in the twist expansion (corresponding to the parton model limit). 
The resulting expression for the asymmetry, known as the Cahn-Gilman (CG) formula [3], 
is given at tree-level by 1 



a rl _ G fQ 9 



/ 20 o \ , ,„ ,1- (1 - y) 2 i 

---^ioL( 1 -T' in8 -') + (i-WM TT ^ f ] • < 9 > 



Here y is the kinematic variable defined as 

IP- (I- £') , . 

where P M , and £' denote the four momenta of the deuteron, the incoming electron, and 
the outgoing electron respectively. In the lab frame, one has y = (E — E')/E where E and 
E' denote of the energies of the incoming and and outgoing electrons. The corrections to 
this Cahn-Gilman formula can be parameterized by writing the asymmetry as 



.4 



HL 



GfQ 2 9 r _ 1 - (1 - y) 



2-, 



ax + a 2 



(11) 



2V27ral0l z l + (l-y) 2 
where the parameters a,j (j = 1,2) are schematically written as 

a,j = -\ (2C ju - C jd ) [l + ^(new) + ^(sea) + ^(CSV) + Rj (TMC) + Rj(B.T)] (12) 

and ^(new), i?j(sea), i?j(CSV), i?j(TMC), and i?j(HT) denote respectively corrections 
arising from possible new physics beyond the SM, sea quark effects, CSV, target mass 
corrections (TMC) [2D arid higher twist (HT) contributions. If one is interested in 
looking for signals of new physics beyond the SM that can leave a footprint in the asymmetry 
via the contributions i?i2(new), it is crucial that all the SM electroweak and hadronic 



corrections to the Cahn-Gilman formula in Eq. (12) are under theoretical and experimental 



control. One can take an alternative viewpoint and instead view a precision measurement of 



Am as a probe of hadronic physics that modifies the Cahn-Gilman formula as in Eqs. (11) 



and (12). 



The analysis of this paper is focused on the higher twist correction Ri(HT) that enters 
the d\ term of the asymmetry. The leading contribution to i?i(HT) appears at twist-four, 
giving rise to a 1/Q 2 power law dependence. In contrast, the leading contribution from 
iJi(TMC), which will also have a 1/Q 2 power law contribution, will be suppressed relative 
to -Ri(HT). The relative suppression of i?i(TMC) can be understood by noting that the 
derivation of the Cahn-Gilman formula is valid even for a finite target mass so that target 



1 We observe that the RHS of Eq.(4.89) of Ref. [H] should be multiplied by (-1). 



6 



mass corrections will always appear in conjunction with at least one of the already small 
effects that correct the Cahn-Gilman formula. 



Given that all the remaining contributions to 5>i in Eq. (12) have at most a logarithmic 
dependence on Q 2 , one can, in principle, make a clean extraction of -Ri(HT) by studing the 
Q 2 dependence of the 5i term in the asymmetry. Similar statements can be made for the 
higher twist effects i?2(HT) that contribute to the a 2 term of the asymmetry. However, the 
study of i?i(HT) is particularly interesting because the leading contribution to i?i(HT) that 
arises at twist-four is given entirely by a single matrix element that characterizes quark-quark 
correlations in the deuteron as first observed in [171 [18]. This is in contrast to ^(HT), which 
receives contributions from several different twist-four matrix elements making it difficult 
to interpret the effect of correlations among quarks and gluons in terms of any one of the 
these matrix elements. 

Before giving the explicit expression for i?i(HT) that we derive below, we first review 
some of the standard notation used in PVDIS pehenomenology. The general expression for 
the asymmetry Am, is given in terms of the five structure functions F^ 2 and F^ 2 3 takes the 
form (21] 

-17 Z ri [1 1 /„, 1 xMl T7<lZ\ 



A 



RL 



( G F Q 2 \ 9 e A (2xyF? z - 2[1 - ljy + ffig) 
Uv^rtJ 2xyF? - 2 [l - 1/y + ] F, 



g e v x{2-y)Ff 



(13) 



This general expression reduces to the Cahn-Gilman formula when the leading twist and 
isospin limits are applied to structure functions and when sea quark and CSV effects are 
ignored. The F^ 2 and F$ z structure functions arise from the interference of the electro- 
magnetic current with the vector and axial part of the weak neutral current respectively. 
The asymmetry is often parameterized in terms of the ratio of the longitudinal to transverse 
virtual neutral vector boson cross-sections 



a 



7 ( 7 Z) 



-2 b 2 



7(7Z) 



7 ( 7 Z) 



a 



2xF 1 



- h 



1 + 



{M 2 x 2 



Q 2 



(14) 



In terms of R 7 ^^ the asymmetry in Eq. (13) takes the form given in Eq. (|2 
quantities Y\ 3 are defined as 



Y-i 



where the 



1 + Rr 



l + (l-y) 2 -y 2 


1 - 


-r 2 /(l + W z ) 


- 2xyM/E 


\ + {l-y) 2 -y 2 


1 


-r 2 /(l + W) 


- 2xyM/E 



Y, 



i-(i-y) s 



l + {l-yf 



l-r 2 /(l + W) -2xyM/E 



(15) 



In this notation, the Y\ and Y 3 terms arise from the interference of the electromagnetic 
current with the vector and axial-vector weak neutal current respectively. 
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One of the main results of this paper is that the relation 

W = W z 



r 2 - 1, 



(16) 



known to hold at leading twist due to the Callan-Gross relations of the structure functions, 
also holds even after the twist-four contributions to i?i(HT) are included at tree level. 
Equivalently, the relation 



(17) 



is valid at twist-four up to perturbative corrections in a s {Q 2 ). However, the the twist-four 
contribution does affect the ratio F? Z /F± in the Yi-term of Eq. fcl) as 



Fx 



1 z 



6 



(2C*i u — Cid) 



CG + HT 



9 



1 



sin 2 # 



l + R^HT) 
+ i?i(HT) 



(18) 



20 

5 V ~ ~9~ 

where the first term corresponds to the Cahn-Gilman limit and where, in the second line we 
have omitted the electroweak radiative corrections for simplicity of presentation as we will 
do throughout much of the remainder of the paper. 



III. ISOLATING THE TWIST-FOUR CONTRIBUTION 
A. Structure Functions 

In this section we review the basic phenomenology and conventions for electron-deuteron 
PVDIS. The differential cross-section for electron-deuteron scattering takes the general form 

d 2 a a 2 E' 



dQdE' 



L 1 

Q4 E \ l*» 7 



GfQ2 u z w% 



(19) 



where E and E' denote the energies of the incoming and outgoing electron respectively in 
the lab frame. The square of the momentum transfer via the exchanged photon or Z-boson 



-[I — f) 2 where and £' denote the four-momenta of the incoming and 



i S Q2 = _ g 2 : 

outgoing electron respectively. The leptonic tensors in Eq. (19) are given by 



Tl 

r-yZ 



(20) 



where A denotes the sign of the initial electron helicity with A = 1,-1 for positive and 



negative helicity states respectively. The hadronic tensors in Eq. (19) take the form 



(x\j;W\D(p)y(x\jj\D(p)) + (x\j;\D(p)r(x\j^\D(p))}, (21) 
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where and denote the quark current coupling to the exchanged photon and Z-boson 
respectively and M denotes the deuteron mass. The hadronic tensors are parameterized in 
terms of the structure functions F^ 2 and ^23 as 

wj. - ( - 1~ + q -f) % + ft - ^V) (p - ^V) Sr q < 

vr* = ( - „„ + «f) % + (p, - (P - ^) Sr g + %Cr F - 

(22) 

The structure functions depend on two variables conventionally taken to be Q 2 and Bjorken 
x = Q 2 / (2P ■ q). The definitions of the F^ z 3 structure functions in terms of the vector and 
axial vector neutral weak current operators can be obtained by first breaking up the weak 
neutral current into its vector (Jy„) and axial-vector (Jf„) parts so that 

Jjx = + J/ln-i (23) 



which allows a decomposition of the hadronic tensor W^ Z in Eq. (21) as 

Wlf = + W$ z , (24) 

where W^g yZ correspond to the hadronic tensors arising from the interference of the 
electromagnetic current with the vector and axial-vector weak neutral current respectively 
and are given by 



K^ nZ = ^ B 27r ) 35(4 W p - ?) 

X 

{X\J Z {A)I ,\D{P)Y{X\J2\D{P)) + (X\j;\D(P))*(X\jZ {A)v \D(P)) 



In terms of the above hadronic tensors the F± z 3 structure functions are given by 
= ( - 9r + ^) + ft - ^) ft " ^) " ^ 



(25) 



q 2 J M V M ? 2 - g 2 ^'MP-q 



^ 2MP-q 3 ' 



(26) 



B. Isospin decomposition of structure functions 



We now show in the limit of good isospin and negligible sea quark contributions, the 
twist-four contributions to the Y 1 term in Eq. (|2| come purely from four-quark twist-four 
operators up to possible higher order perturbative mixing effects involving quark-gluon or 
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purely gluonic operators. This result was first pointed out in [T7J [18] . Here we recast the 
argument in the more modern language, derive an explicit expression for the matrix element 
of the four-quark twist-four operator as a linear combination of the structure functions PJ 7 
and F± z , and provide a corresponding formula for the shift in the asymmetry, Pi(HT). 
Moreover, the matrix elements of four-quark twist-four operators are known [2"H |2"5"] to 
satisfy the Callan-Gross relation. We exploit this property to show that Y\ — 1 up to twist- 
four and that the twist-four contribution in the Y\ term lies entirely in the factor F^ z /F^. 
This result implies that the Y± term in Eq. ^ is in principle a relatively clean probe of 
twist-four quark-quark correlations. 

Following the notation of Ref. [18J, we start with an isospin decomposition of the elec- 
tromagnetic current and the vector part of the WNC as 

J v / = 2[(l-2sin 2 ^)^-^sin 2 ^-(^-^sin 2 MA M ], (27) 
where the isovector, isoscalar, and strange quark currents are, respectively, 

v n = o^V- 4 _ dj^d), = -{u^^u + efyjd), A M = sj^s (28) 

and where we have omitted heavy quark contributions or simplicity (including them is 
straightforward). Using this isospin decomposition of the currents in the expressions for 



W^ v and W^j lZ given in Eqs.(21) and (25) respectively, we arrive at the following isospin 



decomposition for the hadronic tensors 

W^ z = 2(1-2 sin 2 9 w )W;i - ± sin 2 9 W W™ + • ■ ■ , (29) 

where the dots indicate contributions from strange and heavier quarks and the hadronic 
tensors W™ ,ss are defined as 

= ^E( 27r ) 35(4) ^- p -^^i^i^(^))*(^i^i^(^)) ; 



M 

x 

1 



d A x e^( J D(P)|^(a;)^(0)|D(P)), 



2ttM 

\- Y^W^Kpx -P- q)(X\s,\D(P)r(X\s u \D(P)), 



M 

x 



2-kM 



d 4 x e^ x {D{P)\s^x)s v {0)\D{P)). 

(30) 



We ignore subleading contributions arising from the strange and heavier quarks in this 
analysis for simplicity. Contributions to the hadronic tensors involving a product of the 
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isovector current with the isosinglet s v current vanish by isospin symmetry since the 
deuteron is an isoscalar state. 



Next we note that the difference of the and hadronic tensors is given by W t 



[IV 



da 
fiu 



' ' [IV " ' [IV 



w vv 

'' [IV 1 



2nM 



1 , T 



d A x e tq - x (D(P)\-{d(x)-f^d(x) u(0)j v u(0) + (u d)}\D(P)). 



(31) 



As seen above, the operator in is just O (x) of Eq.(J3| which is manifestly a twist-four, 
four-quark operator involving the different up and down flavors of quark bilinears. In the 
context of the light cone operator product expansion (OPE), it contains no local operators 
involving only two quark or two gluon fields as occur at twist-two since the fields located at 
different positions along the light-cone have different flavor. 

The relation in Eq. (31) can be understood from the definitions of W™ and given in 
Eq. (30) and noting that 

Vfj,(x)v v (Q) - s li (x)s„(0) = 



1 rT 



{d(x)^^d(x) m(0)7^m(0) + (u <H- d)}. 



(32) 



We now define flavor- dependent structure functions p^ ss > du corresponding to the hadronic 
tensors W™; ss,du as 



W, 



[IV 



vv,ss,du 



[IV 



r^vv,ss,du t-) t-) t-ivv, ss,du 

(-a +ML)fl + (p -LH Q \(P -^a)^ 



(33) 



so that from Eq. (31) we have the relation 



^1, 



rpss jpdu 

r l,2 ~ r l,2 



(34) 



which allows us to eliminate the F± 2 structure functions in favor of Ff 2 and F d 2 . The 



structure functions Fy 2 can be related to the p^ ss > au structure functions via Eqs. (29) 
and (33) as 



vv,ss,du 



77.7 _ rpw _|_ 2. J?ss 
r l,2 ~ r l,2 * g 1-2 

20 



^1,2 
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rpss rpdu 

r l,2 r l,2J 



2(1 _ ^ S in 2 6 W )F% - 2(1 - 2 sin 2 6 W )F$, 



(35) 



where we have used Eq. (34) to eliminate F^ 2 in favor of F( s 2 and F d 2 . 
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C. Isolating twist-four contribution to the asymmetry 



Using Eq. (j35h for the structure functions Fj^^ that appear in Eq. (|13|), the electron 



polarization asymmetry can be brought into the form 

GpQ 2 2g% (l - f sin 2 9w) F ss - 2g\ (l - 2 sin 2 T du + g e v x(2 - y)F 



A 



RL 



V ! X - (J 



10 



du 



(36) 



where we have introduced the shorthand notation 

F ss = 2xyFl s -2\l - l/y + 

J= du ee 2xyF du - 2 \l - 1/y + 



IT J 
xM 

~E~\ 



^2 > 



rpdu 
b 2 ■ 



(37) 



Next we note that the leading twist contribution to F( s and F^ satisfies the Callan-Gross 
relation so that 

F 2;LT = 2x F?*LTi (38) 

where the subscript LT indicates that this relation generally holds only for the leading twist 
contributions. It has also been shown [24T426] that the four-quark twist-four contribution to 
Ff u and F du satisfies the Callan-Gross relation so that 

F du = 2xF* u . (39) 

We outline an alternate derivation of this Callan-Gross relation for Ff% in Appendix |A~| 



Equations (38) and (39) allow us to write 

•J tt — ^" L,r \-LT 

2xF* u 



LT 
-pdu 



y-2 + 2/y- 
-2 + 2/y 



2xMi 
2xM 



E J 



These relations allow us to write 
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Ft 



du 



~rss 
•'LT 



n. 



(40) 



(41) 



LT 



Using Eqs.(37), (38), (39), and (41) the terms in Eq. (36) proportional to g e A can be brought 
into the form 



A RL 



GpQ 2 2g% (l - f sin 2 9 W ) - 2g A (l - 2 sin 2 9 W ) J 



("4V2 



•ss 
LT 



10, 



9 F du 

10 FfL ■ 



( G F Q 2 \ 9a (i - f sin2 e w) ~ 9% (l ~ 2 sin 2 6 w j 
\2V2na 



du 
s 

LT 



9 Fi 
10 



9 / G F Q 2 



10 V2 v / 27ra 



9a 



20 

1 - — sin 2 6^ ) - 



1 F[ 



du 



10 F( s 



+ ... 



LT 



(42) 
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where we have used the symbol A V RL to denote the part of the asymmetry Am proportional 
to g\ that arises from an interference of the electromagnetic current with the vector weak 



neutral current. The first equality in Eq. (42) is obtained by dividing the numerator and 



denominator of the terms proportional to g\ in Eq. (36) by T ss and using 



du 



da 



-ess 
LT 



+ subleading terms, 



(43) 



to make the replacement =^ — » The subleading terms above denote contributions 

J- LT 

arising from the twist-four matrix element J rdu multiplying subleading twist contributions 



to jF ss . The second equality in Eq. (42) is obtained by using Eq. (41) and the last equality is 



obtained by expanding to linear order in the quantity Ff u /F* S LT . The expression for Ay L in 



Eq. (42) is just the sum of the leading twist Cahn-Gilman term and a twist-four contribution 
from a single four-quark matrix element F\ . Comparing to Eqs. (11) and (12) we obtain 
the main result of this paper: 



fli(HT) 



9 



du 



10 (9 - 20 sin 2 B w ) F[ 



(44) 



LT 



We now derive expressions for F( S LT and Ff u in terms of the F^ structure functions 
which will be useful for phenomenological analyses. The leading twist structure function 
F*? LT can be related to the leading twist F^p structure functions as 

9 _ 1 



F , 



LT 



— F 1 
10 ^ 



2(1 



^sin 2 ^ 



r l;L 



LTi 



(45) 



which follows directly from Eq. (35). At leading twist the structure function F^. LT takes the 
well known form 



F Ilt( x ) = l^2 e2 q ^D(x) + q D (x)), 



(46) 



where qo{x) and q~D(x) denote the quark and anti-quark deuteron PDFs respectively. Treat- 
ing the deuteron as an isoscalar combination of the proton and neutron so that we have the 
relation 



1 



q D {x) 



(q p (x) + q n (x)), 



(47) 



where q p , n (%) denote the proton and neutron PDFs respectively and ignoring sea-quark 
contributions we get 



LT 



X) 



1 

4L 

5 

36 



e 2 u [u p (x) + u n (x)] + e 2 d [d p (x) + d n (x)] 
10 



[u p (x) + d p (x)] 



9 



ss 
LTi 



(48) 
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where we have used the isospin relations u Pjn {x) = d n<p (x) and where the last equality follows 
from Eq. (45). Thus, the leading twist structure function F( S LT can be simply expressed in 



terms of known proton PDFs. Substituting the resulting expression into Eq. (44) leads to 
the result 



fli(HT) 



5(1 -f sin 2 W ) 



da 



Up{x) + dp(x) 



(49) 



Next we turn to the four-quark twist-four contribution Ff u . From Eq. (35), Ff u is given 
in terms of the standard structure functions and F± Z , which can be extracted from 
experiment, as 



du 



(9-20 sin 2 6 w )Fl-hFl 



■>yZ 



(50) 



Note that the LHS of Eq. (50) is manifestly a four-quark twist-four matrix element while the 



RHS includes twist-2 contributions, two-quark twist-four contributions and four-quark twist- 
four contributions. This implies that the twist-2 and two-quark twist-four contributions on 



the RHS cancel out. This allows us to write Eq. (|50|) as 

Tidu 



(9-20sin 2 ^)F 1 7;49 -5F 1 " 



(51) 



where the superscript 4q indicates that only the four-quark twist-four operator contributions 



to F^ are kept. This makes the four-quark twist-four nature of the RHS in Eq. (50) 
manifest. 



Using Eqs.(48) and (51 ) in Eq. (44), the twist-four contribution Ri(HT) to the asymmetry 
from the interference of the electromagnetic current and the vector weak neutral current is 
given by 



fli(HT) 



4 
5 



(9 -20 sin 2 9 W )F1' M - hF[ 



1 - £ sin 2 6 



w 



u p (x) + d p (x) 



(52) 



D. Equality of W and FC< Z at twist-four 



While the results in Eqs. (44) and (49) embody the observation of Refs. [T71 [18] in 



the form of structure functions, the relationship to the parameterization of Eq. ^ for the 
asymmetry is not manifest. To make the implications for the latter apparent, we draw on 
the analysis of the previous section to show that the relation 



(53) 
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is valid at twist-four, implying that Y\ = 1 up to perturbative corrections in a s {Q 2 ). Using 
the following decomposition of the structure functions 



1,2:LT 



+ W: 



(54) 



where F^-lt denotes the leading twist contribution to F^^' and dF^^' denotes the higher 
twist contributions, we can write 



7(72) 



r 



2 r 2;LT ~^ or 2 



2xF$? 
5F i(iZ) 



[1- 



SF, 



7 ( 7 Z) 



1:LT 



-]-! + -, 



7(7^) 



ZXr l:i 



(55) 



where we have expanded F- 7 ^ 7 ^ and the dots denote terms suppressed by higher powers 
of Q 2 . In what follows we only keep terms up to twist-four and suppress + • • • terms in 
Eq.(55). From this relation it follows that 

l + W z 8Ff -2x8Ff 5F2-2X5F? 



l + R-y 



2rF lZ 

z " Lr l;LT 



^ X F\-LT 



(56) 



Using Eq. (35) we can write 
10 



F 1 

r l,2;LT 

5Fl 2 



9 

10 



us 



2;LT' 



F1%LT 2 ( 1 g Sin2 0W)F'i-2:LT- 



8Fr 



rpdu 
P l,2i 



ml 



20 



2(1 - — sin 2 9 W )8F^ 2 - 2(1 - 2 sin 2 9 W )F™,(57) 



where 5F^ S 2 denotes the contribution to F( 2 from terms beyond twist-2. Using the expres- 

(58) 



sions in Eq. (57) in Eq. (56) we arrive at 

{Fi u - 2xFf u ) r 



1 + W Z 



l + Ry 



1 + 



2xF{ 



ss 
LT 



L10 



1 - 2 sin 2 9 
l-^sm 2 9 



w 



w 



Now using Eq. ( 39 ) we arrive at the result 

l + W z 



(59) 



l + Ry 

valid up to twist-four, neglecting perturvative corrections in a s (Q 2 ). This leads to the result 

Y 1 = 1, (60) 

up to corrections in a s (Q 2 ) and power suppressed terms beyond twist-four. Finally, using 
Eqs.(35) and (54), the ratio F± Z / F^ which appears in the Yi term of Eq. (J2J) can be written 

as 



Ff_ 

n 



20 

~9~ 



sin 2 9w) — 



1 F{ 
10 F£ 



du 



LT 



(61) 



leading immediately to Eq. (42). 
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IV. MODEL ESTIMATES OF HIGHER TWIST EFFECTS 

Given the unique sensitivity of the Y\ term in the deuterium asymmetry to the four- 
quark HT operator Ot^(x), it is useful to provide model estimates of its contribution as a 
benchmark for the Jefferson Lab PVDIS program. To that end, we utilize the MIT Bag 
Model[23], following the analysis of Refs. [33] 13%] . In principle, one may consider the use 
of other models to estimate the matrix element of 0^(x), such as QCD sum rules (see, 
e.g. Ref. [35] and references therein) or the instanton vacuum approximation [35] that have 
been applied more extensively to HT effects in polarized structure functions. In addition, a 
nonperturbative QCD computation would yield a result from first principles. To our knowl- 
edge, the particular nucleon matrix element of interest here has not been computed in any 
of these approaches, though some indications may be inferred from related calculations. For 
example, the authors of Ref. [3B] showed that in the instanton vacuum approximation, four 
quark matrix elements are suppressed relative to those of two-quark/gluon HT operators. 
The authors of Ref. [37] have carried out a quenched lattice computation of the contribution 
of the isospin two, four quark operator to the pion structure function using Wilson fermions, 
and find that its scale is set by the square of the pion decay constant, F%. These authors 
suggested that the scale of the corresponding nucleon matrix elements would be set by 
rather than F n , presumably leading to a larger value than in the constituent quark picture. 
In both cases, only contributions to the leading moments were considered. 

In what follows, we will use the MIT Bag Model to estimate not only the overall magnitude 
of the HT four quark contribution but also its dependence on x B . We note that the MIT 
Bag Model was previously employed by the authors of Refs. [201 EH]- m applying their 
computation to the deuterium asymmetry, these authors assumed that the x^-dependence 
of the twist-two and twist-four contributions to the structure functions were similar and 
obtained the twist-four contributions by rescaling the twist-two contributions by the ratio 
of their leading moments. Under these approximations, they obtained 

^^-i^ • < 62 > 

In what follows, we extend the analysis of Ref. |20j by allowing differences in the in- 
dependences of the twist-two and -four contributions to the structure functions. To this 
end, we compute a series of twist-four structure function moments, fit these moments to a 
parameterization in moment number N, and perform an inverse Mellin transform to obtain 
the structure function. This procedure is subject to several uncertainties, including the 
truncation of the tower of moments before carrying out the inverse Mellin transform and 
the choice of parameterization to which they are fit. In addition, we have neglected the 
logarithmic evolution of the moments from the hadronic scale to the Q 2 of interest, as a 
full computation of the anomalous dimension matrix - including the effects mixing between 
the four-quark and twist-four quark-quark gluon and purely gluonic operators remains to be 
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completed. Nonetheless, we believe the computation described below provides a reasonable 
model estimate for the four-quark structure function relevant to i?i(HT). We find that the 
value of the structure function moments decreases rapidly with N, justifying the neglect 
of higher moments in the inversion, and that our results for the inverse Mellin transforms 
do not vary appreciably as we change the parameterization used in fitting them. To the 
extent that the logarithmic Q 2 evolution of the moments at these scales is gentle and that 
the quark-quark correlations embodied in the MIT Bag Model capture the dominant twist- 
four physics, then the estimates described below and illustrated in Fig. [T] should provide a 
reasonable benchmark. 

The four-quark twist-four contributions to i^ 7 ^ are given by[2U [25] 



FT Aq (x B ,Q 2 ) = -Yq^YI / dxd V dz e q e g >U 1 (x,y,z)(A(x,y,z,x B ) + A(y-x,y,y- z,x B ) 

q,q' 

- A(y - x, y, z, x B ) - A(x, y,y-z, x B )) 

+ e q e q ,U 2 {x, y, z) (A(x, y, z, x B ) + A(y -x,y,y-z, x B ) 
+ A(y - x,y,z,x B ) + A(x,y,y - z,x B )) , 

F? Z]iq (x B ,Q 2 ) = — / dxdydz e q g^,Ui(x, y, z) (A(x, y, z, x B ) + A(y - x, y, y - z, x B ) 

q,q' 

- A(y - x, y, z, x B ) - A(x, y,y-z, x B )) 

+ e q g^U 2 (x, y, z) (A(x, y, z, x B ) + A{y -x,y,y-z, x B ) 
+ A(y - x, y, z, x B ) + A(x, y,y-z, x B )) 



(63) 

where e q and gi are, respectively, the quark electric charge and vector coupling to the 
Z-boson with G\ q = —g e A gY/2 and 

(y-x)[z-x) {x-y){z-y) (y - z){x - z) 
and the deuteron four-quark operator matrix elements U\ and U 2 are given by 

U^x,y,z) = -j^ [ — ^— e iA V^-*V^-^D|^(0)^^ 
4A 2 J 2n 2n 2ix 

U 2 (x,y,z) = I ^^^e^e^-^e-^-^(D|^nO)^ 7 5^nH^(H^75t> 9 \An)|D)^ 

4A Z J ZTC ZTC Z7T 

(66) 

where g is the SU(3)c color coupling constant and t a denote the generators of SU(3)c in 
the fundamental representation. We have introduced the hadronic scale A in Eqs. ( |63p5" ) 
in order to express the structure functions in terms of dimensionless functions, though the 
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pi Ki z ) are independent of this scale. The N-th moment of a structure function F(x) is 
defined as 

M(N) = [ dxx N ~ 1 F(x), (67) 
Jo 

and the structure function can be obtained from its N-th moment via the inverse Mellin 
transform 

i rc+ioo 

F(x) = — : / dNx~ N M(N). (68) 

1lX% J c-ioo 

The second, fourth, and sixth moments of the contributions to F^" 1 ^' 49 have been 
computed in the MIT Bag Model in [331 EI]. These results are given in Table, [j] and 
below we give a brief summary of the computation (for more of the moment calculation, see 
Refs. [331 EH)- 

The N-th moment of the structure function is parameterized as 

0» 



M w-Ew- (69) 

i=j 

After the inverse Mellin transform, the corresponding structure function takes the form 



i=2 ^ >' 

The structure functions F^ ' ,4q receive contributions from two terms corresponding to the 
Ui and U2 contributions in Eq.(63) 



prAq _ piAq;Ui _|_ piAq;U 2 

p-/Z;Aq _ p-yZ;iq;Ui _|_ p-iZ;Aq;U 2 

(71) 

For the numerical estimates we use the results of the fit for the U\ and U2 contributions at 
Q 2 = 1 GeV 2 which are given by 

prM,ui „ pyZM-m . ^ = QA5 x ^ flg = _ m x 1Q -4 ^ = 516 J x ^ 

F 7;4 9 ;C/2 _ F 72;4 9; i/2 . ^ = m5 x ^-^ Qg = _ 1197 x ht 4 , a 4 = 1191.4 x 1(T 4 , 

(72) 

In order to determine the correction .R^HT), we substitute these results into the numer- 
ator of Eq. (52). For the denominator we use the CTEQ5 pdfs [39]. The resulting shift 
in the asymmetry, for different representative values of Q 2 , is plotted in Fig. ([!]). For the 
range of xb shown, the correction is largest in the valence region, reaching a magnitude 
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M(N) otF? z ' Aq ' ul (x B ) 


M(2) 
M(4) 
M(6) 




-1.72 x 10" 4 
-0.62 x 1(T 4 



M(N) oiFf' M ' U2 {x) 
M(2) I -34.00 x 10" 4 
M(4) -3.77 x 10~ 4 
M(6) -0.05 x 10" 4 



M{N) ofFZ M ' U1 (x B ) 


M(2) 





M(4) 


-1.71 x 10" 4 


M(6) 


-0.61 x 10" 4 



M(N) of F?' Aq ' U2 {x) 


M(2) 
M(4) 
M(6) 


-33.72 x 10~ 4 
-3.74 x 10" 4 
-0.05 x 10" 4 



TABLE I: The first few moments for the structure functions F^ Z,4q,U1 (x), F^ Z ' 4q,U2 (x), F^ ,iq,U1 (x), 
and F^' 9 ' 2 (x) in the MIT Bag model at Q 2 = 1 GeV 2 . These are updated numbers after correcting 
some numerical errors discovered in |33|. [M] . 



commensurate with that obtained by the authors of Ref. [20] ■ The growth for xb near one 
is a result of the relatively quicker fall off with xb of the pdfs; the twist four contribution in 
the numerator is falling less quickly in this region 2 . One can also perform this analysis with 
alternative parameterizations of the structure function moments. This was done in Ref. [33] 
with similar results and we refer the reader to it for more details. 

For the kinematic range that will be accessible in the planned JLab experiments, the Bag 
Model estimate for the magnitude of the correction i?i(HT) lies below the expected exper- 
imental sensitivity. To the extent that the Y x term can be separated experimentally from 
the remaining contributions to the asymmetry, one could test this Bag Model expectation 
by looking for an appreciable Q 2 dependence in the data. The presence of such a depen- 
dence would point to stronger correlations between quarks of different flavors than implied 
by the Bag Model picture, which correlates the up- and down-quarks largely through the 
confinement radius and the Pauli exclusion principle. On the other hand, the absence of 
large power corrections would imply that the Y\ term can be interpreted primarily in terms 
of the underlying electroweak interactions and/or possible CSV in the parton distributions. 
We comment on the implications for probes of CSV and new physics in the following section. 



V. CHARGE SYMMETRY VIOLATION AND NEW PHYSICS 



To the extent that i?i(HT) is either tiny as suggested by the MIT Bag Model estimates 
or large enough to be extracted utilizing the l/Q 2 -dependence, one may hope to use the 



2 In general, the implementation of the parton model in the threshold region — > 1 introduces theoret- 
ical ambiguities, so we avoid this kinematic regime in computing the relative correction. For a general 
discussion and references, see, e.g., Ref. [10] , 
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FIG. 1: The estimate of R\{HT) as a function of the Bjorken variable x for different values of Q 2 in 
the MIT Bag Model. The curves from the bottom to top correspond to the values Q 2 = 4, 6, 8, 10, 12 
GeV 2 respectively. 

deuterium asymmetry as a probe of CSV and/or new physics. In terms of the former, it has 
recently been suggested that HT contributions to the Y\ term in the deuterium asymmetry 
may be too large and too theoretically uncertain to utilize this term as a probe of CSV [21] ■ 
These suggestions were based on the possibility that W and W z could differ substantially. 
We have shown that at finite Q 2 where a twist expansion is still valid, such a possibility 
cannot apply at twist four since R 1 = W z at this order in the twist expansion, up to 
perturbative corrections. We now compare the MIT Bag Model estimate of -Ri(HT) to 
the CSV correction, ^(CSV). To that end, we follow the parameterization of CSV effects 
utilized in Ref. [21]: 



u p 


= u + 


5u 

y 


dp 


= d + 


5d 

y 


u n 


= d- 


6d 

y 


d n 


= u — 


5u 

y 



(74) 
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In terms of the 5u and 5d one has 



i2i(CSV) 



1 / 2C\ U + C\d\ 3 



2 V 2C\ U — C\d J 10 



8u — Sd 
u + d 



(75) 



The Su and Sd have been constrained by structure function data utilizing the ansatz 



5u — 8d = 2nf(x) 

f(x) = x~ 1/2 (l- x) A (x - 0.0909) , (76) 

with k lying in the range —0.8 < k < +0.65. Detailed phenomenological and theoretical 
analyses of CSV effects can be found in Refs. [2H SU E2] • In Fig- |2j we show the relative 
magnitudes of .Ri(HT) and i?i(CSV) for a representative value of Q 2 = 6 GeV 2 and k 
given by the extremes of the allowed range. We observe that the Bag Model higher twist 
correction is considerably smaller than the possible range for CSV effects. To the extent that 
the Bag Model provides a realistic guide for the magnitude of i?i(HT), a series of precise 
measurements of the leading term in the asymmetry could provide a powerful probe of CSV 
effects. 

The implications for probing new physics via .R^new) are less clear. To be concrete, we 
follow Ref. jl3] and consider new parity-violating contact interactions 

£new = ^ eYl 5 e h v f-ff , (77) 

/ 

where A is the mass scale associated with the new physics, k 2 gives the overall coupling 
strength, and the hy are the specific vector current couplings to each fermion /. Retaining 
only contributions from up- and down-quarks, we obtain the corresponding contribution to 
the correction i?i(new): 

v f-16nn 2 \ /w\2 / 2h v -h v \ 
fll(ll ew)=( — ) ( x ) ( 1 _ 2 ; sm ^ y/9 ) , (78) 

where we have expressed the Fermi constant in terms of the Higgs vaccum expectation value 
v = 246 GeV. 

A given scenario for new physics will determine the specific values of k, A, and the h v . 
For example, E 6 grand unified models contain additional U(l) gauge groups that may lead 
to the existence of a TeV-scale Z' boson. To illustrate the sensitivity of the Vx-term to 
this scenario, we consider a particular pattern of symmetry-breaking that gives rise to a 
low-mass Z x boson. In this case, the correction i?x(new) arises from tree-level exchange of 
the Z x . In terms of the parameters appearing in Eqs. ( 77|78 ) one has k 2 = 2.2a, A = M- 



h v = 0, and h v = -1/20. For M x = 1 TeV, we obtain Rx(new) = 1.85 x 10~ 3 independent 
of Q 2 . Comparing with Figures [T] and [2j we observe that the scale of this correction is 
commensurate with that of i?i(HT) in the MIT Bag Model and well below the allowed 



21 



0.03 



0.02 



0.01 



0.00 



-0.01 



-0.02 



-0.03 



0.3 



k = -0.8 



k = 0.65 



0.4 



0.5 



Ri(CSV) 



Ri(HT) 



Ri(CSV) 



X 



0.6 



0.7 



0.8 



FIG. 2: The relative magnitudes of R\{HT) and R\{CSV) as a function of the Bjorken-x variable 
for a representative value of Q 2 = 6 GeV 2 . using 5u — 5d = 2nf(x) where f{x) = x _1 / 2 (l — x) 4 (x — 
0.0909) for k = —0.8. The top curve and bottom curves give R\{CSV) for the choices k = —0.8 
and k = 0.65 respectively in Eqs.(75) and (76). The middle curve is the MIT Bag Model estimate 
for R^HT). 



bands for the possible CSV correction. In order for a measurement of the l^-term to probe 
this scenario, one would need an experimental sensitivity of ~ 0.2% with knowledge of the 
CSV and HT corrections at a similar or better level theoretical precision. On the other 
hand, the planned 4% measurement of the proton's weak charge by the Q-Weak experiment 
at Jefferson Lab will probe the same scenario for a one TeV Z x . A similar comparison with 
other scenarios suggests that for a determination of i?i(new) to be competitive with the 
Q-Weak experiment as a probe of new physics 3 , one would need a combined experimental 
and theoretical uncertainty of better than ~ 0.5%. At present, then, it appears that a study 



3 We note that the RHS of Eq. (45) of Ref. [43j contains an error and should be multiplied by a factor of 
eight. As a result, the mass bound scale factor for 5\ in Table I should be multiplied by 2^/2. The same 
factor should be applied to the last entries in Tables II-IV. 
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of the Y\ term in the deuterium asymmetry is better suited as a probe of hadron structure 
than of new physics 4 . 

VI. CONCLUSIONS 

Parity-violating electron scattering has become a powerful tool for probing both novel 
aspects of hadronic and nuclear structure as well as possible indirect signatures of physics 
beyond the Standard Model. Its efficacy depends on both significant experimental advances 
in controlling systematic uncertainties and attaining high statistics as well as on substantial 
developments in the theoretical interpretation of the parity-violating asymmetries. PVDIS 
represents a prime example of this synergy between experiment and theory. The first mea- 
surements of the deep inelastic asymmetry for a deuterium target relied on the simplest 
parton-level description of hadrons, yet the result with a 17% experimental uncertainty (for 
the two highest energy points) was sufficient to single out the Standard Model descrip- 
tion of the weak neutral current interaction from other alternatives. Today, one anticipates 
lower-energy measurements at Jefferson Lab with experimental errors below one percent for 
individual kinematic points, making for (9(0.5%) combined uncertainties on quantities of 
interest. The challenge for theory is to provide a framework for interpreting such precise 
results. 

In this study, we have attempted to do so for the leading term in the deuterium asymme- 
try. In principle, it can be kinematically separated from the subleading term (suppressed by 
1 — 4 sin 2 9\y), making it an object of interest in its own right. In going beyond the simplest 
parton model description of the deuteron structure and the Standard Model description of 
the weak neutral current interaction, one may expect contributions to this term arising from 
higher twist operators, the violation of charge symmetry in the leading twist (parton model) 
terms, and new physics. Drawing on early work by Bjorken and Wolfenstein, who showed 
that this term depends on the matrix element of a single, non-local twist four operator, we 
have delineated the twist four contribution from those that may arise from CSV and new 
physics. In doing so, we have shown that that to this order in the twist expansion, one has 
W = W z up to perturbative corrections, making for a theoretically cleaner interpretation 
of the asymmetry than recently suggested in the literature. We have also utilized the MIT 
Bag Model to estimate the x^-dependent twist-four correction and find that it is small com- 
pared to the range of possible CSV effects as implied by global fits to structure function 
data. Typical contributions from new physics are also smaller than the allowed CSV range. 



4 We note that the possible contributions from supersymmeric extensions of the SM have been analyzed 
recently in Ref. [44], though the analysis applied to the asymmetry as a whole and not the Y\ term alone. 
After taking into considerations constraints from other electroweak precision observables and direct search 
limits, corrections of up to 1.5% on the asymmetry are currently allowed in supersymmetric models. 
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To the extent that the Bag Model provides a reasonable guide to higher twist quark-quark 
correlations, one would not expect to observe appreciable sub-leading power dependence on 
Q 2 in the leading term but would, on the other hand, be able to make a clean interpretation 
of this term in terms of CSV. On the other hand, experimental evidence for a substantial 
Q 2 power dependence would point to interesting non-perturbative dynamics underlying the 
higher-twist matrix elements. Either way, a determination of this leading term at the level 
of precision expected for the Jefferson Lab experiments would provide new insights into the 
behavior of non-perturbative QCD. 
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Appendix A: Callan-Gross relation: = 2xFf u 
It has been shown [23H26] that the contribution from the four-quark twist-four operator 



in Eq. (31) to the structure functions satisfies the Callan-Gross relation. Here we 

recast this argument using the language of the Soft-Collinear Effective Theory(SCET) [27- 
[29] . which is an effective field theory for describing the interactions of collinear and soft 
degrees of freedom and can be applied to electron-deuteron scattering in the Breit frame. In 
the language of SCET, the argument for the Callan-Gross relation becomes manifest via the 
structure of the leading order SCET operator that appears at twist-four from a tree-level 
matching. 

Recall that the contribution of the twist-four four-quark operator to the hadronic tensor 

w$ = -^y ^ ^toi^x) mo) + 'OJIW. 



(Al) 



is parameterized in terms of the structure functions Ff% as 



(A2) 



as first written in Eq. ( 33 ) . From this parameterization we note that the Lorentz invariant 
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FIG. 3: Tree level matching of the operator 0^" onto the SCET operator in the Breit frame. The 
dashed fermion lines indicate collinear fields in the standard notation of SCET. 



quantity P^P v Wf^ is given by 



in' 



p u wf du 



q 2 M 



da 



T?du -, 

2x 



(A3) 



If it can be shown that this quantity vanishes, it implies the Callan-Gross relation 

2xF? u . (A4) 



We formulate the argument in the Breit frame where the momentum of the virtual photon 
or Z-boson is 

W — 



<r 



and the momentum of the deuteron is 



Q- 



» m 2 

p» = fi ■ P— + . 

2 fi-P 2 



(A5) 



(A6) 



where we have introduced the light-cone vectors 



1,0,0,1), n" = (1,0,0,-1), n 2 = n 2 = 0, n ■ n = 2. (A7) 
From these relations and using momentum conservation in the Breit frame we have 

(A8) 



n-P M 2 

~ <C 1 

n-P~ Q 2 ' 



so that the deuteron momentum is entirely along the light-cone up to power corrections 
in M 2 /Q 2 . 

In SCET in the Breit frame, at leading order in the power counting in Aq CD /Q 2 and at 
tree level, the four-quark operator in Wj™ will be matched onto an SCET operator [15] 

d r y fl d(x) u^ v u{Q) = n^n u / duJidiO^diiJ^duj^ C(ui, u; 2 , u; 3 , u^) 



(A9) 
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where we have used standard SCET notation, a detailed explanation of which, can be 
found in [271429] . This is schematically shown in Fig. [3j Here C(ui, o>2, <^3, ^4) denotes the 
matching Wilson coefficient. W denotes the momentum space Wilson lines 



W 



perms 



cxp 



(A10) 



and the fields A njP , £„ iP denote collinear gluon and quark fields which are fourier transformed 
to momentum space with respect to large light cone momentum component in the n M di- 
rection. The Wilson lines W, determined by collinear gauge invariance in SCET, generate 
all the spin terms at twist-four. The labels p, Ui on the SCET fields denote the large part 
of the light cone momentum and the x dependence of the collinear gluon and quark fields 
corresponds to residual momentum fluctuations. For more details of the matching and the 
SCET notation used here see [45J. What is relevant to our discussion is that the collinear 
quark fields £ ra that appear in the SCET operator satisfy the equation of motion 



4in = 0, 



(All) 



which is the reason that there are no terms proportional to n fl n l/ in the tree level matching 
in Eq. (A9). Thus, contracting Eq. (A9) with n^n v gives zero using the property n 2 = 0. 



Applying these considerations to Eq. (Al) we get the relation 



(A12) 



which from Eq. (|A3j) implies the tree level Callan-Gross relation of Eq. (A4) up to power 
corrections in Aq CD /Q 2 . 
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